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The flavour-changing neutral current transition 6 — ► s can serve as a sensitive probe 
of WIMP dark matter models, if the WIMP mass is under 2 GeV. In this work we extend 
our earlier analysis to a generic class of models where the interaction between the dark 
matter sector and the Standard Model matter sector is mediated by the Higgs boson(s). 
We show that experimental limits on the decays of S-mesons to K{K*) and missing 
energy provide stringent constraints on the parameter space of such models, but do not 
rule out sub-GeV WIMPs in a model-independent way. We find that in the context of 
the NMSSM with light pseudoscalar Higgs, the WIMP masses under a few hundred MeV 
are generically excluded with the exception of few highly tuned points in the parameter 
space. 
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1. Introduction 

The dark matter content of the total energy density in the Universe is now well- 
measured, with several experiments determining an abundance of rtoAih^ ^ 0.12 
So far the presence of dark matter was deduced only through its gravitational 
interaction, and thus its identity remains a mystery. 

Among numerous candidates for dark matter suggested in the literarture, the 
weakly interacting massive particles (WIMPs) are the most interesting from the 
point of view of particle physics. The residual abundance of such particles in the 
early Universe is well understood, with the observed dark matter abundance pro- 
viding a measure of the WIMP couplings through their annihilation cross section at 
freeze-out (for a review, see Ref. Furthermore the inferred couplings of WIMPs 
to the Standard Model fields are sufficiently large that detection of WIMP pair- 
production via a missing energy signal is possible in current particle physics exper- 
iments. 
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The mass of the WIMP can be constrained from above and below using certain 
theoretical arguments. The assumptions that the mass scale of the particles that 
mediate interactions between Standard Model and WIMPs is not lighter than the 
weak scale, and couplings are perturbative produce the WIMP annihilation cross 
sections that vanish both in the limits of very large and very small WIMP masses. 
Comparison with (JannVrei ~ 0.1 pb , which is required to reduce the dark mat- 
ter abundance to an acceptable level during freeze-out, produces the allowed mass 
range, referred to as the Lee- Weinberg window ^"^ Obviously, such a window is 
very model dependent. In many popular models this limit requires the WIMP mass 
to be larger than a few GeV. However it is possible to construct simple models in 
which the required masses can be as low as 100 MeV (e.g. Ref 1181 1191 120(1 . The 
Lee- Weinberg limit is usually quoted for fermions, in which case the annihilation 
cross-section is suppressed by a factor of mj^j^j/M'^, where M is the mass of the 
mediator particle, which requires large couplings for light WIMPs. However if the 
cross section is enhanced, or if scalars are used instead of fermions, the lower limit 
on WIMP masses can be significantly reduced. 

Light dark matter models have several interesting features. The dedicated dark 
matter searches, such as DAMA and CDMS, rely on measuring the recoil of nuclei 
from WIMP scattering and as such are less sensitive to WIMPs which are signifi- 
cantly lighter than the nuclei used. Light WIMPs produce small recoil energy, which 
is usually well under the background that peaks at small E. As indicated in Figure^ 
this insensitivity to sub-GeV WIMPs will not be solved with the next generation of 
underground experiments, although the CRESST experiment is expected to provide 
an improved probe of dark matter as light as 0(1 GeV) due to its low recoil thresh- 
old. Interestingly enough, there is a possibility that a light WIMP could explain the 
positive signal observed in the annual modulation at DAMA without violating the 
bounds from other experiments 

It has also been suggested that the annihilation of light WIMPs could explain the 
galactic positron excess*^. The positron excess has been measured t hrou gh the flux 
of 511 KeV 7-rays meassured by the INTEGRAL/SPI experiment l^^^, and seems 
to require non-standard sources of positrons. Although there are other possible 
astrophysical explanations (see for example Ref. |H1 IHl ll()|l , positron production by 
annihilations or decays of light dark matter particles 1^ remains as one of the more 
attractive solutions due to its natural explanation of the nonlocalized nature of the 
observed 7-rays. It has also been demonstrated that the annihilation of ligh t dark 
matter does not conflict with measurements of either the galactic 7-ray fluxEEl 
or the extragalactic 7-ray backgroun , and may also provide an explanation for 
the observed flux of 7-rays in the 1-20 MeV range^. 

As we recently demonstrated, an intriguing possibility of a sub-GeV WIMPs 
is the opportunity of detecting a pair of WIMPs as a missing energy signal in B 
meson decays^! Indeed, as dedicated underground searches for dark matter remain 
largely insensitive to sub-GeV WIMPs, their production is possible in experiments 
such as BaBar and BELLE, which produce large numbers of i3-mesons, and can 
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Fig. 1. Current limits on WIMP-nucleon spin-independent cross sections from dedicated under- 
ground searches. The solid lines represent the predictions for the minimal scalar model with a 
100 GeV Higgs, while the current limits are given from (I) CRESST, (II) CDMS (2005 Si), and 
(III) CDMS (2005 Ge). In the interval of 100 MeV - 2 GeV the predicted signal has signiciant 
QCD-relatcd uncertainty. 

study their rare decay modes. As a result such facilities provide a new opportunity to 
search for light dark matter. For the minimal scalar WIMP model these experiments 
have already excluded most of the parameter space with ms < 1 GeV, while future 
data from B factories will be able to probe as high as mg ~ 2 GeV - ' . 

The purpose of this paper is to explore the question of how generic the limits on 
light WIMPs derived in Ref. El^re, and whether all dark matter models with sub- 
GeV WIMPs can be efficiently constrained by i3-physics. To answer these questions 
we study the class of models where the interaction between Standard Model sector 
and WIMPs is mediated by one or more Higgs particles. We demonstrate that 6 — > s 
decays with missing energy provide important constraints on the parameter space 
of such models. We also point out the possibility, based on the two-Higgs doublet 
model (2IIDM) at large tan/3, that these constraints can be circumvented. 

In Section 13 we review our previous results on the minimal scalar model and 
extend the result for more general scalar models with an additional singlet scalar 
that mixes with the Higgs boson. In Section 01 we apply the same tecniques to a 
related model with two Higgs doublets and calculate the branching ratios of WIMP- 
producing decays of B-mesons. This model has the additional benefit of relaxing the 
fine tuning condition required for a sub-GeV scalar WIMP in the minimal model. In 
Section^ we introduce some simple models of fermionic dark matter, calculate the 
WIMP production in i3-decays, and discuss the limitations on such models from 
the Lee- Weinberg limit. We also address the case of NMMSM (next-to-minimal 
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supersymmetric Standard Model) where we show that WIMPs under 200 MeV 
are generically excluded except special points in the parameter space with 6 — > s 
transitions being fine-tuned well below its natural scale. In Section we review 
the experimental bounds on the decay B ^ K + missing energy, and discuss the 
potential for future searches. 

Before we present our calculations, we would like to make several important 
comments. In our discussions of the Lee- Weinberg limit, it should be noted that 
there is no clear division between the perturbative and non-perturbative regions 
of parameter space and therefore no definitive lower bound on the WIMP mass. 
However in this paper we will require the WIMPs to be sufficiently heavy that the 
abundance constraints satisfy < 47r, with k representing an effective coupling in 
each model. 

As it is done in Ref. 1171 we treat the significant hadronic uncertainties in the 
annihilation cross section using the most "optimistic" and "pessimistic" scenarios. 
That results in the certain range for predicted WIMP-producing branching ratios. 
In most of the models we present, the largest uncertainty exists for a WIMP with 
mass Too m ~ 500 MeV that could annihilate through the fo resonance potentially 
enhancing the cross section and resulting in weaker couplings. However, it is unclear 
whether such resonances exist at the freeze-out temperature of T ~ 0(50 MeV). The 
heavier resonances are not expected to contribute to the annihilation cross section 
as the higher freeze-out temperatures will significantly reduce the strength of the 
resonances. We also will assume that WIMPs heaver than ~ 1 GeV will freeze-out 
before the transition from quarks to hadrons. 

We will also limit our discussion to the decay channel B K + missing energy. 
While it is possible to produce WIMPs in other decays, such as i? — > missing energy 
or B —^ J + missing energy, the experimental limits on such decays are inferior 
to B ~* K + missing energy. Heavier WIMPs could be produced in T ^ 7 + 
missing energy, however a SM decay of T is due to strong interactions, which 
greatly reduces the branching ratio of WIMP production compared to weak B- 
meson decays. 

Finally, it should be noted that the collection of models presented in this paper is 
not exhaustive. There are other possibilities which could result in light dark matter 
[19 20 21^ that generically require relatively light particle mediating the interaction 
between the dark matter sector and the SM sector. In this paper, we will breifly 
touch a possibility of a lighter pseudoscalar Higgs particle discussed in Ref. |2] and 
omit the models with an extra U{1) gauge group that have an additional Z' in 
the sub-GeV range. The latter possibility is interesting, and should such Z' have 
axial-vector couplings to fermions, one expects large flavour-changing b ^ s — Z' 
amplitudes induced by radiative corrections. However, should this new force couple 
to a pure vector current (with e.g. B — L quantum numbers), the constriants from 
flavour-changing decays will be greatly relaxed. The detailed analysis of models with 
light WIMPs supplemented by light Z' goes outside the scope of the present paper. 
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2. Minimal Scalar Models 

The simplest WIMP model is a sing let scalaiE2E3IHI 

which interacts with the Stan- 
dard Model through exchange of the Higgs: 



(1) 

= X^' + + + ^VEwS'^h + ^S''h\ 

where H is the SM Higgs field doublet, vew — 246 GeV is the Higgs vacuum expec- 
tation value (vev) and h is the corresponding physical Higgs, H = (0, (vEw+h) / \/2). 
The physical mass of the scalar S receives contributions from two terms, rrig = 
uiq + Xv%]Y, and requires significant fine-tuning to provide a sub-GeV mass. In this 
section we will calculate the branching ratio for the pair production of scalars in the 
decay B ^ K + SS, which contributes to Br(_B+ + missing energy). Being 

minimal, this model obviously possesses maximum predictivity, and the branching 
ratio of WIMP production can be calculated as a function of dark matter mass only. 

It should be noted that the decay B ^ K -\- missing energy is actually ex- 
pected to occur regardless of the existence or nature of light dark matter. As shown 
in Figure Et- and |2b, the Standard Model predicts the transition h ^ s + Vv 
at one loop, so that the B-meson can decay to neutrinos with Br(i?+ 

+ missing energy) ~ (4 ± 1) x 10^^. However as demonstrated before^! 
the decay B ^ K + SS (resulting from the b s transition shown in Figure ^) 
can enhance the missing energy signal by up to two orders of magnitude. 

The transition b s + h occurs as a loop process, which at low momentum 
transfer can be calculated by differentiation of the 6 — > s self-energy operator with 
respect to vew, 



( 'ig^mbmlV*^Vtb \ _ 
^bsh = ^,0,^2 ...... ^LbBh + (h.c.). (2) 



As the Higgs is significantly heavier than the other particles involved in the process, 
it can be integrated out leaving an effective Lagrangian for the & — > s transitions 
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with missing energy: 



(3) 



Leading order Wilson coefficients for the transitions with dark matter scalars or 
neutrinos in the final state are given by 



C. 



DM 



m 



h 



327r2 



(4) 



9w 9w^t*sVtb 



167r2 



C( + 2xt ixl — 6x 



8{xt ~ f) 8(xt - 1)2 



In xt 



where Xt = m'l/M'^. From this effec tive th eory and the hadronic form factors 
calculated in the light-cone sum rules l ^^ l '^^ l ^ the partial width for decays with 
missing energy can be calculated 



Br, 



*K+SS 



~ 4 X f 0" 
where we use the parametrization 



2.8 X lQ-'^K^F{ms). 



\ mh J 

and the available phase space as a function of the unknown tos, 



(5) 



(6) 



F{ms) 



h{sfl{s,ms) ds 



h{sYl{s,Q) ds 



with 



I{s,ms) = [f - 2s(M% + Ml) + {Ml - Mlf]^\l ~ 4^2/5] i 

For comparison, BaBar recently reportecP^ a limit of Br(i?+ + 
missing energy) < 5.2 x f 0^^ at 90% c.l. which is well above the SM prediction. 
Belle's prelimin ary results were reported at Br(i3"'" — > K'^ + missing energy) < 
3.5 X 10"^ level 1251. Similar calculations can be used for the decay B K*SS, 



BiB+^K+' + K- 1-3 X 10"^ + 3.0 X IQ-'^K^Fims). 



(7) 



with an analogous form factor. 

Even before we perform the freczc-out abundance calculation and extract k, it 
is obvious that the decay into WIMPs would dominate the Standard Model rate 
to neutrinos by one-to-two orders of magnitude near the Lee- Weinberg bound of 
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(a) Constraints on k from dark matter abundance. The parameter space above A 
produces a density too small to explain dark matter, while the space below B would 
lead to a higher than observed dark matter abundance, (b) Predicted branching 
rat ios for the decay B'^ K~^-\- missing energy^ with current limits from CLEO 
(I) EH , BaBar (II) EHl, Belle (III) and expected results from BaBar (IV). 
Parameter space above curves I, II and III is excluded. The grey bar shows the 
expected B Kvv signal. Parameter space to the left of the vertical dashed line 
can also be probed with K'^ tt"'" + missing energy. 
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K ~ 0(1). To make this statement more precise, we extract the parameter k using 
the relation of the WIMP annihilation cross section, 



'annVrel = ^^^X 1™ f^^) (8) 



with the observed dark matter abundanceD, VlBMh"^ ~ 0.12. Here Th-,x is the total 
decay width for a virtual Higgs of mass '2ms, and in our calculation we utilize 

the zero temper ature deca y widths which were studied previ ous ly in early searches 
for a light Higgs l ^ l '^'^ l ^^ l as explained in our previous paper The allowed range 
of K is plotted in Figure[5^, where the main uncertainty comes from Th^hadrons- It 
should be kept in mind that the domain of k's &eZow curve B gives an over-production 
of dark matter and therefore is firmly excluded. Inserting k — ms domain allowed 
by the relic abundance into the result for the branching ratio gives the prediction 
of Br(_B+ — > + missing energy) as a function of ms, FigurelSt. 

These results can be easily generalized to models of scalar dark matter coupled 
to the SM via an additional singlet Higgs particle U. The simplest model of this 
type has the potential 



-Cs = ^S' + + ii^iU + ^i2U^)S^ + V{U) + i^'U^H^H 

(9) 

2 2 

where in the second line we retained only mass terms and relevant interaction terms. 
u denotes the excitation around the vev of U and fj, and rj stand for dimensionful 
parameters presumably of order the electroweak scale. The last term in the second 
line of Q gives the mixing between scalars u and h. If such mixing is significant, 
the existing bounds on the higgs mass would also place a lower bound on the mass 
of the u-boson. The effective Lagrangian for b s + f! transitions is the same as 
Eq. with the Wilson coefficient given by 



Cdm - — — 2 — — \^^) 



which, with the redefinition 



2 ^ /^ lOO GeV y 

~ ml \ mu J 

results in the same abundance constraints and the same branching ratio for B 
K + missing energy, as plotted in Figure |21 Thus, we see that the model with 
more parameters in the singlet sector @ in the limit of light WIMPs gives identical 
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predictions to the minimal model and in this sense the predictions of Figure 
are generic. 

3. Scalar dark matter in 2HDM Model 

The singlet scalar model of WIMPs coupled to the SM is the most economical model 
of dark matter, and for sub-GeV WIMPs can be well constrained with existing B- 
factory experiments. However, the light masses of WIMPs in this model are not 
natural, as it requires significant fine-tuned cancelation between ttiq and Xv'^yy. 
However, if the electroweak sector of the SM is modified, this fine-tuning can be 
significantly relaxed. 

In this section, we consider a singlet scalar WIMP which interacts with two 
Higgs doublets, Hu and Hd- 

_£^!^52^Ai52(|i/°P + |i7,-p) + A,5^(|i?°P + |i/+P) + A35^(iJ,-i/,+ -i7M 

(11) 

We assume the most conservative flavour arrangement in which the up- and down- 
types of quarks originated from the expectation values of different Higgses, exactly 
as it happens in supersymmetric models (see e.g. Ref. I30|l. This model differs from 
the SM in a significant way: due to two different vev's of the two Higgs doublets, 
there is an additional parameter tan/3 = Vu/vd, which is large when Yukawa cou- 
plings in the down-sector are enhanced. Indeed, there are many theories which 
attempt to unify the Yukawa couplings of the third generation of the Standard 
Model (eg see Ref. |33 ESI ■ Such unification typically requires tan /3 w mt/mt or 
tan/3 « mt/mr and thus Vd ~ 0(few GeV). Therefore in the regime of large tan/3 
the mass corrections which depend only on Vd will in general be small and may lead 
to minimal fine tuning. 

For example, the physical mass of the WIMP scalar in the model given above is 

rrig = ml + Xivj + X2vl - XsVuVd. 

If tan/3 ~ 0(1), sub-GeV WIMPs would require the same amount of fine tuning as 
before. However in the special case of large tan /3 and the hierarchy of A's, Ai ^ 
A2,A3, the mass correction depends only on Vd, resulting in ~ 0(1 — 10) for 
a GeV-scale WIMP mass, which then does not require significant fine tuning. In 
the case of A3 » Ai, A2, the fine-tuning is also relaxed, especially if the abundance 
constraint allows for A3 0(tan~^/3). In the rest of this section, we analyze the 
WIMP production in the model, assuming large tan /3. 

3.1. Ai dominant 

The additional diagrams which contribute to the decay B K + SS and give the 
leading-order contributions in tan/3 are given in Figure |31 The effective Lagrangian 
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Fig. 3. Diagrams contributing to the decay 6 — > s + SS in the 2HDM plus scalar dark matter 
model when Ai is dominant and tan /? is large. Inside the loops, and denote the two charged 
Higgs bosons, with the mixing of the two doublets denoted by a cross. 



is of the same form as Eq O and in the limit of large tan (3 and ^ the 
Wilson coefficient is 



r, ^1 gw^tsVtbXt ( I - at + adnat ^ . 

^"''"W. 32.2 (1_,^)2 ) (12) 



where at — mj/Mfj. The branching ratio is calculated as before, 



Btb^k+e - 4.0 X f 0-« + 3.2 x IQ-'^^^ ( ^ (t^"^)^"* ) ^i^s) (13) 



2 



Here we use the parameterization, 



Ai 



100 GeV\ 
Mh J 



2 



and unlike the minimal case there is an additional rather mild dependence on jtih 
through at- ft disappears in the limit mH S> mt. We also make a safe assumption 
that additional diagrams with charged Higgses are not going to alter the branching 
to neutrinos. As before, F{ms) is constructed such that 

F{ms) = [ '"""^ /o(s)2/(s, ms)ds f '"""^ fo{sfl{s, 0)ds 

and F{0) = 1, and F{ms) = for ms > ^{mB—mK). The branching ratio is plotted 
in Figurelit with the current experimental bounds from BaBarI^, BELLeI^ and 
CLEO-^. 

The abundance constraint on k is calculated as in Sectional with one important 
difference. When tan (3 is large, the scalars predominantly couple to leptons and 
down-type quarks. Consequently, the charm and top loops do not contribute in 
the effective couplings to giuons, leading to the reduction in the cross-sections for 
SS <?5,7r7r by a factor of 0(l/iV^) ^ 0.1, where Nh is the number of heavy 
quark flavors that convert virtual Higgs to hadrons. Nh = 3 in the minimal model, 
and Nh = 1 in the 2HDM with large tan/3. The charm threshold also does not offer 
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Fig. 4. (a) Constraints on k from the observed dark matter abundance (b) Branching Ratios for 
B ^ K -{- missing energy in the two higgs doublet model, with WIMPs coupled primarily to /fj. 
The labeling of current limits from BaBar, BELLE, and CLEO is the same as as in Figure 1^. 



additional annihilation channels. As a result, the values of Ai that would fit the 
observed abundance of dark matter are higher than the corresponding values for A 
in the minimal scalar model studied previously. The resulting k, Figure clearly 
corresponds to a strong interaction regime below WIMP masses of 500 MeV. 
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Fig. 5. Diagrams contributing to the decay b ^ s + SS in the 2HDM plus scalar dark matter 
model when A3 is dominant and tan /3 is large. 



3.2. A2 dominant 

The case in which the scalars couple predominantly to produces results similar 
to the minimal model of Section [2| In the large tan/3 limit Vu ~ vsm and therefore 
the fine-tuning of the scalar mass is significant, while the branching ratio for B —^ 
K + SS and constraints on A2 are not expected to be significantly different from 
the minimal model with one Higgs doublet. 



3.3. A3 dominant 

The third possibility is that the dark matter scalars couple to the HuHd combination 
of the two Higgs fields. In this case the mass of the scalar is 

2 2 ^3,V%m 
Too = TOn — 

^ " tan/3 

which in the large tan/3 limit requires less fine-tuning to produce sub-GeV WIMPs 
than in the minimal scalar model. 

We find that in this model the branching ratio of WIMP pair production in 
B-meson decays is suppressed. It turns out that two b ^ s + SS diagrams that have 
tan (3 enhancement shown in Figure |31 exactly cancel each other so that Cdm ~ 
O(tan0 13). 

At the same time, the abundance constraint is similiar to the Ai dominant case, 
except that the scalar annihilation cross-section is enhanced by a factor of tan^ /3, 
resulting in A3 being tan (3 times smaller than Ai from the previous example. A com- 
bination of the enhanced annihilation cross section and cancellation of ^ s + SS 
amplitude at leading order in tan/3, results in tan~^ (3 ^ ©(lO^"^ — 10^'*) suppression 
of 6 — > s + SS branching ratio relative to previous examples. Consequently, 



BiB^K+ss < BiB^K+uiy for A3 dominant (14) 

and therefore this model cannot be constrained by i?-decays. It is then clear that 
a generic case with both couplings, Ai and A3, being important can create a large 
range for J^tb^k+ss even for a fixed value of tos, and therefore the model of scalar 
dark matter coupled to 2HDM is far less predictive than the minimal model. 
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Fig. 6. The Feynman diagrams which contribute to b ^ s + f! in the 2HDM plus ferinionic 
WIMP model. 

It also interesting to note that the elastic scattering of dark matter on nucleons 
is enhanced by tan^ /3 when A3 is dominant, and therefore this model will indeed 
create a measurable signal for detectors with low recoil energy threshold such as 
CRESST. Detailed analysis of such constraints goes outside the scope of the present 
paper. 



4. Fermionic Dark Matter 

In the previous two sections the WIMPs were presumed to be scalar fields. This is 
common in models of light dark matter, as the Lee- Weinberg limit on the WIMP 
mass is in general lower for scalar WIMPs. However fermionic dark matter could 
still be light if its annihilation cross section is enhanced. In this section we will 
present some simple models of fermionic dark matter, and review the constraints 
from _B-decays. 

The simplest renormalizable model of fermionic dark matter is the analogue of 
the scalar model given in EqlHl 



(15) 



ml .,2 



-YXX + + + VVewwUh + XuU^ 



where x is a Majorana fermion and H is the Standard Model Higgs field. As before, 
w =< U >~ 0{vew)- However the annihilation cross-section, at freeze out, in this 
model is suppressed relative to the scalar model by 



^scalar 



0(10-5) 



^There is a possibility of producing light fermionic WIMPs in this model if ruu <S feio, and 
r) <Si li < 1, though such a model would require significant fine-tuning to avoid being detected in 
direct searches of light(er) Higgs. 
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Fig. 7. Abundance constraints on k in the 2HDM plus fermionic WIMP model. It is expected 
that most of the parameter space near curve B is also excluded, however a detailed calculation 
of the annihilation cross-section is beyond the scope of this paper, and as such we have 

utilized the perturbative formula for production of gluons and unbound quarks. However except for 
possible resonances the decay is expected to be significantly smaller, and therefore the abundance 
constraints would require the coupling constant to be larger. 



for ~ 0(1 GeV) and 111^ ^ 0{vew)- As a result of this suppression, the fermions 
cannot annihilate efficiently in the early Universe and therefore the couplings would 
have to be non-perturbative to explain dark matter. 

The second model we will consider is the analogue of the model presented in 
Section|3| For the model to be renormalizable, the Majorana fermions cannot couple 
directly to the higgs fields, but must instead couple through an intermediate scalar. 

+ \2U'i\Hy + \H+\') + hU'{H^H+ H'.Hl) + X^U' 

As in the previous model, the inclusion of fermions in this model requires an increase 
in the coupling constants by a factor of ~ m'^/{wm^), with w —< U >, to produce 
the same abundance of dark matter. However in the special case of A3 S> Ai, A2 the 
annihilation cross section is also enhanced by a factor of tan^ and so x could be 
light without requiring large A3. The abundance constraints are given in Figure [3 
with 



2 .x2 2 f^smwY f 100 GeV Y /tan 
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Fig. 8. Feynman diagram for the process b — > s + XX iii the final model. 

As in the A3 dominant case presented in Section |31 the 0(tan/3) contributions (see 
Figure 0) to the process b ^ s -\- xx cancel and as a result 

Br(B ^ K + xx)) < Br(S ^ isT + I7i/). 

The final model we will consider is a Majorana fermion x coupled to higgs- 
higgsino pairs. By "higgsinos" we mean an SU{2) x U{1) charged fermionic fields 
with the same quantum numbers as the Higgs fields but do not impose supersym- 
metry requirements on the size of the couplings. By the same token x can be called 
the neutralino. The relevant part of the Lagrangian is given by 



where M <C /.t, Xu^m and tan/3 is large. In this model the dark matter candidate is 
taken to be the lightest mass eigenstate, 

X = -Tpcose + HdSinO ain^ 9 = 



(17) 



mi = M 1 - 



It should be noted that this model is constrained by the Z-boson invisible decay 
width, and requires sin^ < 0.15. 

The b s transition in this model is calculated from a single diagram, given in 
Figure |S1 The effective Lagrangian is 



^b^sf = '^CoMrribSLbRXX - C^SLltibL^ltiV + [h.c.) (18) 



with Wilson coefficient 



_ V't*Vtfctan/3 / ArfA„w„^ A a* In a* 
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Fig. 9. (a) Constraints on k from the observed dark matter abundance. {b)Branching Ratios for 
B ^ K + missing energy in the neutralino modeL Current Hmits from BaBar, BELLE, and 
CLEO are indicated as in Figure I^Jj. The Standard Model predictions are given in the gray box. 



where at = mj/Mfj, and the branching ratio is 



BRB^K+f:^4,-0xl0-^ + 9.8x10-^ K^^^] F{mi), (20) 



2 



where 
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100 GeV 



y / tan/3 
j 100 



The phase space integral for fermions is defined as 



fo{sYis-2m'/M'B)I{s,ni^)ds 



/o(s)'s/(J,0) drs 



(21) 



The cosmological abundance and the branching ratio are plotted in Figure |2| with 
the current experimental limits. Although it appears that the effects of dark matter 
in this model could be observed, the abundance constraints require k to be in the 
non-perturbative region for < 1 GeV. Furthermore, constraints on A„ from Z- 
decays and on Mh^ require the original coupling constant Xd to be a few orders of 
magntude larger than k and therefore the model is non-perturbative for < 2 
GeV. Heavier WIMPs may still be produced in B-decays, though the branching 
ratio is significantly smaller than the Standard Model signal. 



5. NMSSM with light dark matter 

One of the most popular dark matter candidates is the lightest neutralino present 
in supersymmetric models. The existence of supersymmetry is well motivated in 
particle physics, and the cosmological abundance of neutralinos is comparable to 
the dark matter abundance for a large region of parameter space. In this section we 
will demonstrate that light neutralinos, if kinematically allowed, can be produced 
in _B-meson decays and will provide a significant contribution to the branching ratio 
Br{B ^ K + f]). 

In discussions of supersymmetric dark matter, it is common to work with the 
minimal supersymmetric standard model (MSSM) . However in that case the lightest 
neutralino is expected to be heavier than a few tens of GeV and therefore too heavy 
to produce in _B-meson decays. Although it is possible to lower this limit to 6 
GeVl^, it requires extreme tuning of the parameters, while sub-GeV WIMPs seem 
to be impossible even with the fine tuning. However it is possible to produce sub- 
GeV neutralinos in extensions of MSSM by introducing new resonances significantly 
lighter than Mw- 

In this section we present a case study of NMSSM with light neutralino dark 
matter. The lightness of x in this niodel is made possible, without significant fine- 
tuningl^, by one of the pseudoscalar Higgs bosons a being chosen to have a mass of 
a few GeV. A comprehensive analysis of dark matter in this model was attempted 
in Ref. |^ while the constraints from _B-physics were considered in Ref. |^ None 
of these papers, however, considers the decays of _B-mesons with missing energy. 
As mentioned in the previous section, the neutralino annihilation cross section in 
the non-perturbative QCD regime has not been adequately addressed. A proper 
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Fig. 10. Additional Feynman diagram for the process b —> s + XX the NMSSM model. 

account of < 1 GeV would require the calculation oi a ^ gg ^ PPP and 
a ^ qq ^ PPP, where P is a pseudoscalar meson such as tt and K, plus an 
extensive scan over the NMSSM parameter space, which we will not attempt here. 
Instead, we choose to analyze the domain of neutralino masses of < < 1.5m7r 
where the annihilation proceeds primarily via niuons in the final state, and is free 
from hadronic uncertainties: 

2 / X 2 



^ ann ^re 



1 / Aq^ sin 7 tan /3 \ / mj 
2tt \ ml J \VEW 



-rei - ( ->^-r"^ ] [7^] ™y 1 " -^/-^- (22) 



In this formula, the coupling constants Xax a-nd mixing angles 7 originate from the 
tree-level interaction Lagrangian of a with the SM fermions and neutralinos, 

- Qnt = T^^ax^xilsX + siu 7 tan /3 ^^0/175 /. (23) 
Z ^-^ Vew 

Here the summation goes over the charged leptons and down- type quarks, 7 
parametrizes the admixture of IvciHj^ in a, and tan [3 is assumed to be a large pa- 
rameter. Specific relations between 7, the neutralino-pseudscalar coupling Aq^ and 
the parameters of the fundamental NMSSM Lagrangian can be worked out, but are 
of no interest in this discussion since they bear no consequences for our predictions. 

The SUSY diagrams that can play an important role in the 6 ^ s -I- xX de- 
cay are shown in Figure 1101 These are two Higgsino-stop exchange diagrams that 
have different parametric tan/3 and 7 dependence. They give the folowing Wilson 
coefficient in the effective Lagrangian defined in (|18|l . 

r, Vj* l/tbg^xt tan /3 /Aqx sin 7 tan /3\ 

'-'DM = ITZ-^ 2 ^SUSY- (i4J 

In this expression, we adhere to the minimal flavour-violation scheme of the soft- 
breaking sector and treat the stop mixing as a mass insertion. The loop function 
I SUSY depends on the trilinear soft-breaking parameter At, masses of left- and 
right-handed stops, which for simplicity can be chosen equal, rrij^j^ ~ '^Ijit — ; 
and /i parameter: 

Atfi f 1 — af + aj:hi 
IsusY = — 2" n Y2 " 



1,1 -^aff 



tan /3 sin 7 



(25) 



where 
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/i^/m? Note that in the context of NMSSM the Higgsino-pseudoscalar 



couphng Ajjjj and /i-parameter are often chosen to to be related via the singlet's 
vev. It is important that Isusy does not decouple even in the limit of large SUSY 
masses, provided that the relative size of ^t, fi and rrif remains fixed. When all 
SUSY masses are equal, At — fi = rrif, the expression in front of the square bracket 
in H25|l becomes 1/2, which we would consider as a natural value for Isusy- Notice 
an additional power of tan/3 in the expression for Cdm I|24|I . a very well-known 
enhancement in the SUSY models ^'"^ . By introducing an effective coupling k, 



tanpY flO GeV" 



10 y V m, 
we express the i3-decay misisng energy signal as 



(Xax sin 7)^ 



BiB^K+f = 4.0 X 10-s + 0.98^2 ( 



SUSY 



with 



F{m^)^ / h{sfsl{s,m^)ds 



H -1 



h{sfsI{s,Q) ds 



(26) 



(27) 



(28) 



Using the freeze-out abundance to constrain k. 



< 



2.4 GeV 



(29) 



and imposing the experimental bound of Br(i3+ + < 3.5 x 10^^ on the 

branching ration, we arrive at the following constraint on Isusy, 



/ tan (3 



Isusy] < 3.0 x 10""* for = 150 MeV 



(30) 



This is a requirement on Isusy to be suppressed relative to its natural value by 
three orders of magnitude. We conclude that unless an unnatural cancellation of 
Isusy occurs, the possibility of a ~ 100 MeV neutralinos in NMSSM is excluded. 
However, even if the SUSY contribution vanishes for some reasons, the next-to- 
leading order in tan (3 contribution mediated by charged Higgs loops is still capable 
of producing an observable signal, which from existing bounds still excludes a ~ 100 
MeV neutralino. 



6. Experiment 

The decays of the type + f! discussed here involve missing energy but 

provide no constraint on the associated missing mass. This is problematic experi- 
mentally, even for a hypothetical detector with full solid angle coverage, unless one 
of the two B mesons in the event is fully reconstructed. Similarly, it is hard to 
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imagine that a sensitive search can be done without the well-constrainted initial 
state provided by experiments at the threshold for BB production. This discussion 
will therefore focus on experiments at the T(4S'). 

The B reconstruction can be done using either hadronic or semileptonic decays. 
After the decay products of one B are removed from consideration, signal events 
should contain only a single kaon while background events will generally contain 
additional charged or neutral particles. In reality, beam-related backgrounds tend 
to leave additional energy in the electromagnetic calorimeter and interactions of 
particles from the fully reconstructed B with material in the detector can lead to 
additional, spurious charged or neutral particles. The selection criteria balance the 
loss of efficiency from strict vetoes on additional activity against the suppression of 
background. After removing particles associated with the fully reconstructed B, no 
charged tracks are allowed beyond the identified kaon, and the neutral electromag- 
netic energy deposition is restricted to bo loss than 200-300 McV. The background 
is negligible for K energies near tob/2 and rises steadily toward lower K momenta. 
As a result, the experiments place a requirement on the minimum accepted K mo- 
mentum; this is further discussed below. 

The CLEO measurement reconstructs B mesons in hadronic decays and requires 
Pk > 0.7 GeV, providing sensitivity to dark matter particle masses up to 2.1 GeV 
and placing a limit of Br(B+ K+i^V) < 24 x lO"'"' for the SM decay. The BaBar 
and BELLE measurements each require pK > 1.2 GeV, thereby restricting sensitiv- 
ity to dark matter particle masses below 1.9 GeV. BaBar reconstructs B mesons in 
both hadronic and semileptonic modes with comparable sensitivity in each, and de- 
termines Br{B+ K+uV) < 5.2 x 10"^ based on a sample of 88 x 10^ T(45) BB 
decays. BELLE uses only events with one B reconstructed in a hadronic decay mode 
and find the preliminary result Br(B+ — > K'^vV) < 3.5 x 10~^ based on a sample 
of 275 X 10*^ T(4S) -> BB decays. 

The expected background from normal B decays with two or more missing 
particles is several times higher than the signal expected from the SM process with 
neutrino pair. This background has contributions from semileptonic decays (e.g. 
B Dip with D K~^£~i' or D ^ K'^X) where the the detectable particles 
other than the lie outside the detector acceptance, as well as from hadronic 
B decays involving one or more K^, which have a non-negligible probability of 
traversing the electromagnetic calorimeter without depositing significant energy. 
The latter decays in some cases still have poorly known branching ratios, but these 
can be measured (or more stringently bounded) with larger S-factory data sets. 
The sensitivity can be expected to improve as 1/VN as more data are analyzed. 

The decay B^ + f! doesn't add much sensitivity due to the smaller re- 

construction efficiency for both the (about 1/2 that of _B+) and the /-T" (only 
1/3 decay to tt+tt^). The decays B K* fl have sensitivity comparable to 
— > /r+ -|- f (the branching ratios for both the SM and New Physics transitions 
are higher due to the higher probability of producing a K* in the fragmentation of 
the s quark, offsetting the lower experimental efficiency). The B K* -|-.^ channel 
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should be sensitive to dark matter masses only slightly smaller than those explored 
in the _B+ + f channel, allowing for an improvement in overall sensitivity 

to light dark matter. This, coupled with the ~ 1 ab""'^ data samples that will soon 
be available from the B factories, should push the sensitivity to missing energy 
processes to 2 times the level of the SM branching ratio. 



7. Conclusions 

In conclusion, we have demonstrated that B-decays with missing energy provide a 
viable method for exploring light dark matter models. This alternative to dedicated 
underground experiments has considerable sensitivity to light WIMPs, and can 
probe dark matter using existing experimental results. While this result has been 
previously demonstrated in Refll7| this paper extends the result to include several 
different models. 

In particular it has been demonstrated that the 2IIDM+scalar model, in the case 
of Ai dominant, can satisfy the abundance constraint with Br(i?+ + SS) > 

Br(i?+ + vv). This model has an additional benefit, as it requires only 

minimal fine tuning to achieve sub-GeV WIMPs in the large tan (3 limit. Data 
from BaBar and BELLE excludes the parameter space below ms ^ 1.3 GeV, 
while future measurements are expected to probe up to ms ~ 1.9 GeV. The case 
of A3 dominant also provides sub-GeV WIMPs with minimal fine-tuning. However 
in that case the branching ratio Br(_B+ + SS) is suppressed by factor of 

tan^ /3 relative to the Ai dominant case, and therefore is significantly smaller than 
the Standard Model signal. In this case, only the direct searches with the low recoil 
energy threshold has a chance of detecting these WIMPs. 

Furthermore we have demonstrated that it is possible to search for fermionic dark 
matter in B-decays, though in the models presented the couplings are required to 
be non-perturbative to meet the abundance constraint for < 2 GeV. While it is 
possible that heavier fermions could be produced in B-decays, the branching ratio 
for such a process is significantly smaller than the Standard Model predictions. 
It is also possible to lower the Lee- Weinberg limit in these models by allowing 
other light particles which are very weakly coupled to the Standard Model, such 
as in the NMSSM with one light pseudoscalar Higgs. In such a model, it would be 
possible to have weak couplings and still produce an observable increase in Br(_B 
K -\- missing energy). Our analysis shows that indeed O(IOO) MeV neutralinos of 
NMSSM are generically excluded by the B decay data, while heavier (above 1 GeV) 
neutralinos would probably be compatible with existing constraints. 

Except for the case of NMSSM, we have limited the discussion to models which 
do not include new forces below the electroweak scale. Models with new resonances 
below Mw can be explored using _B-decays, but in addition they might be well con- 
strained by existing experiments and therefore should be introduced with caution. 
Such models can be analyzed in exactly the same manner as presented in this paper. 
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